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A possible mechanism of electronic phase separation in the systems with orbital ordering is an- 
alyzed. We suggest a simple model taking into account an interplay between the delocalization of 
charge carriers introduced by doping and the cooperative ordering of local lattice distortions. The 
proposed mechanism is quite similar to the double exchange usually invoked for interpretation of 
phase separation in doped magnetic oxides like manganites, but can be efficient even in the absence 
of any magnetic ordering. It is demonstrated that the delocalized charge carriers favor the formation 
of nanoscale inhomogeneities with the orbital structure different from that in the undoped material. 
The directional character of orbitals leads to inhomogeneities of different shapes and sizes. 
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I. INTRODUCTION 



The existence of superstructures is a characteristic fea- 
ture of magnetic oxides, in particular those containing 
ions with orbital degeneracy, i.e., Jahn- Teller (JT) ions. 
In the crystal lattice, the JT ions usually give rise to the 
orbital ordering (00)i^. The 00 is typical of insulating 
compounds. The electron or hole doping can destroy 00 
since the itinerant charge carriers favor the formation of a 
metallic state without 00. However, at low doping level, 
we have a competition between the charge localization 
and metallicity. It is well known that such a competition 
can lead to the so-called electronic phase separation (PS) 
with nanoscale inhomogeneities^i^. This phenomenon is 
often observed, e.g., in doped manganites and is usually 
related to some specific type of magnetic ordering, an- 
tifcrromagnetic insulator versus ferromagnetic metal. In 
the usual treatment of PS, the 00 is not taken into ac- 
count (see, however, the discussion concerning isolated 
orbital and magnetic polarons^'^'*'^). Here we study the 
effect of 00 on PS employing minimal models includ- 
ing itinerant charge carriers at the 00 background, and 
show that at small doping the PS may appear in systems 
with orbital degeneracy even without taking into account 
magnetic structure. We consider this effect using two 
versions of the models. First, in Section II we study a 
symmetrical model analogous to the Kondo-lattice model 
in the double exchange limit, where the orbital variables 
play a role of local spins. Namely, it is supposed that 
localized electrons create lattice distortions, leading to 
the formation of 00. The conduction electrons or holes, 
introduced by doping, move on 00 background. In the 
second version (Section III), we take into account the 
specific symmetry of eg type for doped electrons. For 
both versions, we demonstrate the possible instability of 
a homogeneous ground state against the formation of in- 
homogeneities. As a result, additional charge carriers 



introduced by doping favor the formation of nanoscale 
inhomogeneities with the orbital structure different from 
that in the undoped material. In Section IV, we deter- 
mine the shapes and sizes of such inhomogeneities and 
demonstrate that depending on the ratio of the electron 
hopping integral t and the interorbital coupling energy J, 
the shape can vary drastically. For the two-dimensional 
case, in particular, there exists a critical value oft/ J, cor- 
responding to the abrupt transition from nearly circular 
to needle-like inhomogeneities. This is a specific feature 
of orbital case: the directional character of orbitals brings 
about the unusual and very rich characteristics of inho- 
mogeneous states. 



II. SYMMETRICAL MODEL 

Let us consider the system with JT ions having double- 
degenerate state. This degeneracy can be lifted by lo- 
cal lattice distortions, giving rise to two different ground 
states of each ion, a or b (e.g. a (b) state corresponds to 
elongation (compression) of anion octahedra) . The states 
a and b of the ion n determine the corresponding orbital 
states of a charge carrier at this ion. In general case, each 
ion can be characterized by a linear combination of basis 
a and b states, described by an angle 9 



\e) = cos-|a) 



sin-|6) . 
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The local distortions can interact with each other leading 
to some regular structure. In the simplest symmetrical 
case, the interaction Hamiltonian can be written in a 
Heisenberg-like form 
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where Tn = {t^, t^} are the Pauli matrices, and a and b 
states of the ion n correspond to eigenvectors of operators 
r^, with eigenvalues 1 and —1, respectively. For Hamil- 
tonian two simplest kinds of ordering are possible: 
ferro-00 (the same state at each site) and antiferro-00 
(alternating states at neighboring sites). In the absence 
of the charge carriers, the ground state is antiferro-00 if 
J > and ferro-00 if J < 0. Of course, in real materials 
with Jahn- Teller ions, the orbital Hamiltonians are more 
complicated, but the analysis based on the model ^ 
seems to be sufficient to reproduce the essential physics 
related to the orbital ordering. 

Under doping, itinerant charge carriers appear in the 
system, so the density of charge carries n ^ 0. We assume 
that the charge carriers are doped into double-degenerate 
states and move on the 00 background determined by 
localized electrons. (As we argue below, the main results 
will be also applicable to the case where the same elec- 
trons, e.g. eg electrons, are responsible both for the 00 
and for conduction due to doping into these eg states). 
The values of electron hopping integrals should depend 
on the states of the neighboring lattice sites. The electron 
Hamiltonian can be written as 

Hcl = ~ {PLaalaaamf3aPmf3a + h.C.) , 

(nm) .OL,(3.a 

(3) 

where, a\^aa^ o.naa are creation and annihilation opera- 
tors for the charge carriers at site n with spin projection 
(T at orbital a. Having in mind that we are dealing with a 
strongly correlated electron system, we introduced in ([3]) 
projection operators P excluding a double occupation of 
lattice sites (we consider the case n < \). Below, ana- 
lyzing the electron contribution to the total energy, we 
shall consider square lattice in the two-dimensional (2D) 
case and cubic lattice in the three-dimensional (3D) case 
using the tight-binding approximation. For the spectrum 
of charge carriers, we have 

1 ^ 

£;"^(k) = -t'^'^-^Yl = t"^C(k) , (4) 

i=l 

where D is the space dimensionality and ki are the com- 
ponents of wave vector k in the units of inverse lattice 
constant 1/d. 

In our model, doped electrons at a JT distorted site 
a or 6 are in the corresponding orbital state \a) or \b). 
We can introduce three hopping integrals: f", t''^, and 
^ab _ ^ba _ ^1 -p^j. giniphcity, let us assume that t"''^ — 
t^^ — t. Then, we have a competition of two factors: 
the formation either of a wider electron band or of an 
optimum 00 type. 

At the site n in the state 9, the charge carrier has an 
orbital state 16* >, described by Eq. ([1]). The hopping 
integral between the sites characterized by orbital states 
l^i) and 16*2) can be written as 



First, we consider a homogeneous state, assuming that 
the orbital structure corresponds to the alternation of 
l^i) and 1^2) orbitals. In the mean-field approximation, 
we can represent the total energy per site as 

Etot{eue2) = zt'^'^^eoin) + ^ cos(0i - 62) , ea[n) < , 

(6) 

where z is the number of nearest neighbors, and dimen- 
sionless kinetic energy eo{n) is determined by the type of 
the crystal lattice. A specific form of eo{n) for different 
cases will be discussed below. We assume in this Section 
that t, t', and, therefore, t^^^^, do not depend on the 
direction of hopping. In this isotropic case, So{n) does 
not depend on 9i^2, and we can easily calculate the or- 
bital structure by minimization of total energy ([6]) with 
respect to angles 9i and ^2- At relatively large dop- 
ing, when t|£o(?^)| > 2 J, we have ferro-00 state with 
di — O2 — n/2. In the opposite case, i|eo("^)| < 2J, the 
minimization yields: 62 = tt — 0i, and 

<,..arcsi„(^feM), i^<l, (7) 

The total energy of such a canted orbitally ordered state 
is 

Etot^zt'eoin)-^el{n)^^. (8) 

Note that if eo{n) — nf{n), where f{n) varies slowly 
with n, then Etot can have a negative curvature, at least 
at small n, which is a signature of an instability of a 
homogeneous orbitally ordered state (negative compress- 
ibility) . 

Let us now determine function £0(7^) and analyze the 
dependence of the total energy on doping. For the tight- 
binding spectrum ^ of electrons in the lattice of the 
dimension D, the density of states po{E), has the form 

oc 

= / (If^^''"^^''^^ ^ cos(£;.)Jo"(^), 



(9) 

where Jq is the Bessel function. Then we have 

eo{n) = J ^ dEEpoiE), (10) 

with the chemical potential fi given by equation n — 
j!^^dEpo{E). 

At small doping, n ^ 1, it is possible to write eo{n) in 
a simple explicit form. In 2D case, eo{n) « —n + 7rn^/2, 
and the total energy then reads 

, f t'^ Trt'\ y zJ 

Et.t^-z^n-z[^---y--. (11) 
From Eq. ([TT]), we find that d^Etotjdr? < if 
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This implies an instability of the homogeneous orbitally 
canted state toward the phase separation into phases with 
ferro- and antiferro-orbital ordering. The situation here 
is quite similar to that for the usual double exchang&iS, 
which corresponds to the case t' = 0. At relatively large 
t', when 2-Kt' jt > t/J, a homogeneous state is stable in 
the whole range of doping. 

Taking SQ{n) « —n in Eq. ([7]) we get a rough estimate 
for a region of phase separation: 



< 71 < 



2J 



(13) 



So, the orbitally canted state turns out to be unsta- 
ble nearly in the whole range of n where the difference 
02 - Oi = TT - 201 with 01 from Eq. Q is non-zero. The 
situation remains qualitatively the same, if in Eq. ([8]) 
for Etot we take So(n) calculated using the density of 
states ([9]). The behavior of Etot{n) in 2D case is illus- 
trated in Fig. [T] 

In three dimensions, the situation is more complicated. 
At small doping, we have £o("-) ~ + an^/'^ ^ where 
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and the total energy becomes 



Figure 1: (Color online) Two types of behavior of the energy 
of homogeneous state ((8]) in 2D as function of doping n: with 
a region of negative curvature (red solid line, J /t = 0.05), and 
without it (blue dashed line, J/t = 0.35); t' — 0.5t for both 
curves. For the values of parameters corresponding to the 
solid curve, the orbitally canted state existing on the left-hand 
side of vertical line is unstable toward a phase separation. The 
Maxwell construction in a region of phase separation is shown 
by green dot-dashed line. The homogeneous state correspond- 
ing to the blue dashed curve is stable in the whole range of 
doping. 
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The second derivative of Etot is positive at n 
changes sign at 



5aJt' 



9^2 



(15) 



Taking into account the same arguments as in 2D case, 
we get an estimate for the phase-separation range 



<n < 



2J 
T' 



(16) 



We see that the presence of nonzero nondiagonal hopping 
t' leads to the appearance of a lower critical concentration 
Uc for phase separation. (Maxwell construction would 
lead to phase separation is a somewhat broader doping 
range, starting from some hq < Uc). 

Note that inequalities (|13p and are valid at rela- 
tively small values of J/t ratio. 



III. ANISOTROPIC MODEL 



Now we study a more realistic model of eg orbitals 
on the square 2D lattice. This situation is character- 
istic, for example, for layered cuprates, like K2CUF4, 
or manganites (La2Mn04 or La2Mn207). We assume 
that an orbital exchange Hamiltonian has Heisenberg- 



be written as a linear combination of two basis func- 



tions a >= 



> and \b >- 



|2z2 



0, but it \e >= cos(6'/2)|a;2 - > + sin(6'/2)|2z2 - 



y2 >: 
t/2 >. 

The hopping integrals t"^ in Eq. ([3]) now depend on the 
direction of hopping, and can be written in the form of 
a matrix 



3 



(17) 



where minus (plus) sign corresponds to x (y) direction of 
hopping. 

Assuming again an underlying orbital structure corre- 
sponding to the alternation of \0i) and 16*2) orbitals, we 
obtain the spectrum of charge carriers in the form 



E{k) = -to iAxi0i,92) cosk^ + Ay{0i,e2) cosk 



where 



Ax,y{&l,d2) 



(18) 



1 - ^2 \ , V1+U2 .TT 

COSl^j-HCOS^i- 



The total energy then reads 



(19) 



Etot{0U02) = to{A^{01,02)+ Ay{0x,B2))y^ 

ein;0i,02) + 2Jcos{0i- 02), (20) 



like form In the case of eg orbitals, any orbital can where e{n; 0i,02) — J_i dE Ep{n; 0i, 6*2), and the density 
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of states can be written as 



p(n;e 



oc 

/ds , . 
— cos(i?s) X 
TT 



(21) 
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Note that the density of states now depends on angles 
6*1, 62 via functions Ax y{9i,92). In order to find orbital 
structure, one should minimize -Etot, Eq. pp]) . with re- 
spect to 9i and 02- The analysis shows, that at doping 
n less than some critical value ni , depending on the ra- 
tio J /to, the minimum of the total energy corresponds 
to 6*1 = 0, 6*2 = TT, that is, we have the homogeneous 
antifcrro-orbital structure with alternating |a 



2/2 > 



and \2z^ 



' y > orbitals (we ignore here anharmonic 



effects and higher-order interactions, which usually sta- 
bilize locally elongated octahedra with the angles, in our 
notation, 6 = n, ±27r/3, see Refs. [Tllfl^ . The energy of 
such a state is 



(22) 



This state is locally stable, d'^Etot/dv? > 0. 

At n = rii, a jump-like transition to the canted state 
with 02 = —01 occurs, where 



t()\e(){n)\ 
4J 



(23) 



and Etot{n) has a kink at n 
canted state at n > ni is 



ni. The energy of such 



Etot^ta£Q{n)~ ^el{n)-2J . (24) 

With the further growth of n, the angle 9i decreases, and 
ain — n2, determined by the equation io|eo("-2)|/4J — 1, 
it vanishes, 9i = (ferro-00 with \x'^ — y^) orbitals). 
The total energy of the system as function of doping is 
shown in Fig. [51 Note, that depending on the values of 
parameters, the energy (|24p can have either positive or 
negative curvature (see the inset to Fig. [2]). In the for- 
mer case, the homogeneous state is locally stable in the 
whole range of doping, but the phase separation still ex- 
ists in the range of n near n = ni, due to the kink in the 
system energy. In the second case, PS, of course, also 
exists (we have an instability in some range of doping, 
where d'^ Etot / dn^ < 0). Note, that these two possible 
situations (negative curvature of Etot and the kink) can 
lead to inhomogeneous states with quite different prop- 
ertiesi^. 



IV. INHOMOGENEITIES IN THE ORBITALLY 
ORDERED STRUCTURES 

We demonstrated above that the additional charge car- 
riers introduced to the orbitally ordered structures can 
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Figure 2; (Color online) The energy of homogeneous state for 
the anisotropic model as a functions of doping at J /to = 0.025 
(red solid curve). In the region near ni « 0.08, the homoge- 
neous state is unstable toward a phase separation. In the 
inset, the dependence of Etot{n) + toVSn (linear term of the 
dependence Etot{n) in the range n < ni is subtracted) on 
doping n is shown at the range near ni at different model 
parameters. The red solid curve (blue dashed curve) corre- 
sponds to J /to — 0.01 {J /to = 0.035), and have a negative 
(positive) value of d^Etot/dn^ in the region n > ni close to 
111 . The phase separation exists for both situations. Maxwell 
construction is shown by dot-dashed line. 



lead to the formation of an inhomogeneous state. Now, 
let us discuss possible types of such inhomogeneities in 
more detail using a model of the eg orbitals at the sites 
of 2D square lattice, considered in the previous Section. 
We assume that each charge carrier forms a finite region 
of an 00 structure with alternating 16*1) and 16*2) orbitals 
(not necessarily ferro-00 with 9i = 92 = 0) to optimize 
Hc\- The remaining part of the crystal has antiferro-00 
structure with \x'^ — y'^ > and \2z'^ — — y^ > orbitals, 
according to the results of the previous Section at n — > 0. 

The spectrum of charge carriers is given by Eq. (|18p . 
Expanding this spectrum in power series of k up to 
the second order, we find an effective Hamiltonian for 
a charge carrier in a finite region: 



iief! = -to{A.^^Ay) + 



A. 



92 



A, 



92 



^ 9a;2 ^ dy 



(25) 



where A^,, Ay are given by Eq. p9)) . Using Hamilto- 



nian (|25p , we can solve the Schrodinger equation within a 
finite region, which we choose in the shape of ellipse with 
semiaxes \J A^po and ^ AyPQ. As a result, we find the 
following expression for the kinetic energy of the charge 
carrier within such droplet 



Ekin — —to {Ax + Ay) + 



to j 0,1 



(26) 
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to po, we find 




0.5 

(a). J /to = 0.005 




(b) J/to = 0.02 

Figure 3: Total energy as function of angles 9i and 62 at J/to 
smaller (a) and larger (b) than the critical value ~ 0.0075. 



where jo,i — 2.405 is the first root of Bessel function Jq. 
The potential energy Epot related to the orbital order- 
ing is the sum of two contributions proportional to the 
droplet volume v {v = ix ^jA^^Ayp'^): the energy of the 
canted 00 within the droplet is zJcos(0i — Q-i)!'^ and 
the loss in energy of the antiferro-00 matrix due to the 
formation of the droplet is zJ/2. As a result, we get 
(z = 4) 



2 / Cl — C2 



(27) 



Po 



1/4 



(28) 



The total energy (per lattice site) then reads 

Etot = -2J + E{9i, 62)71, (29) 
E{eu92) = ~2taiA,,+Ay) (30) 



+jo,i (8^TtoJ^/A^A^ 



1/2 



where we assume that all charge carriers introduced by 
doping form such identical 00 droplets. 

To find possible types of 00 droplets, we minimize 
E with respect to 9i and 62 (note again, that functions 
Ax,y depend on 61, 62 according to Eq. (fT9|) ). The func- 
tion E{9i, O2) at two different values of J/io is shown in 
FigO In general case, the function E{9\, 62) has several 
minima, and the values of 0i and 82 corresponding to the 
lowest minimum depend drastically on parameter J/tg. 
At small J/to (Fig [3^), the lowest minimum corresponds 
to 9i = 02 = 0, that is, we have ferro-00 structure in- 
side the droplet with occupied |a;^ — > orbitals. In 
this case, the most favorable shape of droplets is a circle 
(see left panel of Fig. [5]). At J/to larger than some crit- 
ical value ( Jcr/to — 0.0075), the minimum 9i — 62 = 
becomes metastable, and the energy E{9i, 62) has four 
degenerate lowest minima; two of them correspond to 
6»i = 7r/3, 02 = 27r/3, 61 = -7r/3, 6*2 = -27r/3, and 
the similar two minima with the replacement 9i 62 
(see Fig[3}D). In this case, we have chains of alternat- 



mg \x — 
and \2x'^ 



> and \2y^ 



> (or \y^ 



y 



>) orbitals and hence nearly one- 



Minimizing the droplet energy Ekin + Epot with respect 



dimensional (cigar-shaped) droplets stretched along y or 
X axes (right panel of Fig. |4]). With the further growth 
of J/to the metastable state 9i = 62 = splits into two 
states corresponding to = —62 with positive and neg- 
ative as it can be seen from Fig The droplets 
corresponding to these states have circular shape, but 
canted 00 structure. 

The existence of two types of droplets with different 
shapes can be easily understood. The maximum gain 
in the kinetic energy corresponds to the ferro-00 state 
with \x'^ — y'^ > orbitals. At small J/to, the kinetic energy 
prevails, and we have circular droplets with this type of 
orbitals. The minimum cost in the potential energy cor- 
responds to nearly one-dimensional structures. At larger 
J/to, the potential energy plays more important role than 
the kinetic one, and we get cigar-like droplets (smaller 
volume of such a droplet gives smaller loss of orbital en- 
ergy). The orbital structure inside the droplet described 
above corresponds to the maximum gain in the kinetic 
energy for one-dimensional chain (in the absence of hop- 
ping between neighboring chains). 

The analysis shows, that the energy of an inhomoge- 
neous state, Eq. P5|) . consisting of circular or cigar- like 
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Figure 4: (Color online) Schematic illustration of circular (left panel) and needle-like (right panel) droplets. An electron or 
hole moves in a finite region creating a ferro- or canted-OO structure, in antiferro-OO matrix. In the case of hole, there exists 
a one mobile empty site within a droplet. 



00 droplets embedded into an antiferro-OO matrix is 
less than the energy of a homogeneous state in a certain 
range of doping < n < n*. With the growth of the 
number of charge carriers, the droplets start to overlap, 
and at n = n* the inhomogeneous state of considered 
type (ferro-00 droplets in antiferro-OO matrix) disap- 
pears. However, the phase separation exists in a wider 
range of doping (see the previous Section) . For circular 
droplets, we have an estimate n* ^ I/tt/Oq. Taking for 
estimate the ratio J/t^ — 0.005, we get po ~ 2 (in units 
of lattice constant) and n* ^ 0.08. 

In the case of cigar-like droplets, we have — (or 
Ay = 0), and according to Eq. ((28|) we would get that 
chains have infinite length (but zero volume v), po — 
oo. This is of course not a very realistic result, coming 
from an approximation, where the potential is assumed 
to be proportional to the droplet volume only. In order 
to estimate the characteristic length L of the chain, we 
should take into account the surface term (proportional 
to the droplet 's length) in the potential energy Epot of 
the droplet. Let us consider, for definiteness, the chain 
of — > and \2x^ — y"^ — z"^ > orbitals, stretched 
along X axis. In this case, we have A^ — VS, Ay — 0. 
The effective Hamiltonian (|25|) is reduced to 



characteristic length of the chain: 



2 dx^ 



and the kinetic energy of the charge carrier in the chain 
of length L becomes Ekm = —toV3{l — 7r^/2L^). The 
surface energy of interorbital exchange interaction has 
a minimum, when the chain is located in an antiferro- 
OO matrix like shown in Fig. [D each — > 
(|2a;^ — — z"^ >) orbital in the chain has its nearest 
neighbor \x'^ — y"^ > (|2z^ — x^ — y^ >) orbital in the ma- 
trix. In continuum approximation, the potential energy 
can be written as Epot = dJL/A. Minimizing Ekin + Epot 
with respect to L, we arrive at the following formula for 



1/3 



9J 



(31) 



At J/to = 0.05 > Jcr/^Oi we have Lq « 5.5. At ran- 
dom distribution of the chains in the matrix (we have 
chains stretched both along x and y axes), the critical 
concentration is about n* ~ ^/Lq, but it can be larger 
if a more complicated structure of chains, e.g. regular 
stripes, appears in the system. 



V. CONCLUSIONS 

We have studied a simple model of electronic phase 
separation in the system of charge carriers moving in an 
orbitally ordered background. It was shown that a ho- 
mogeneous state in such a system can be unstable toward 
a phase separation, where delocalized charge carriers fa- 
vor the formation of nanoscale inhomogeneities with the 
orbital structure different from that in the undoped mate- 
rial. The shapes and sizes of such inhomogeneities were 
determined for 2D lattice of eg orbitals. The shape of 
inhomogeneities depends drastically on the ratio of in- 
terorbital exchange interaction and a hopping amplitude 
of the charge carriers, J/to: there exists a critical value 
of J/to, corresponding to the transition from the circu- 
lar inhomogeneities to a one-dimensional chains of finite 
length. 

The model under study is quite similar to the dou- 
ble exchange model, where the orbital variables play a 
role of local spins. It is well known that such a model 
also exhibits an instability toward a phase separation into 
phases with different types of magnetic ordering. The 
inhomogeneous state with circular ferro-00 droplets is, 
in essence, an analog of a magnetic polaron state (ferro- 
magnetic droplets in an antiferromagnetic matrix) , which 
is usually considered in the double exchange model^i^ii^. 



7 



Nevertheless, our orbital model is more complicated than 
the usual double exchange due to the existence of non- 
diagonal hopping amplitudes and to the anisotropy in 
hoppings. Both these features lead to the results specific 
for the orbital model, such as the kink in the energy of a 
homogeneous state and canted-00 needle-like droplets. 

In the present paper, any magnetic structure and spins 
of the charge carriers were fully neglected. Taking into 
account spin degrees of freedom can lead to the forma- 
tion of inhomogeneities with different orbital and spin 
configurations. 

In the proposed model, the localized electrons forming 
an orbital order and the conduction electrons or holes 
were supposed to be two different groups of electrons. 
However, we can argue that our main results are also 
valid for a model, where the same electrons take part 
both in the hopping and in the formation of orbital or- 
dered structure. Indeed, in the case of magnetic ox- 
ides with Jahn- Teller ions, an orbital degeneracy is lifted 
by lattice distortions, giving rise to an orbitally-ordered 
ground state at n = 1. If we suppose that a long- 
range orbital ordering still exists at small hole doping 
a:: = l — n<Cl, we come to the situation considered 
in present paper: we have holes moving in an orbitally- 
ordered background. In a mean-field approximation, we 
only should replace in all formulas above n ^ x — 1 — n 



and J — > J{1 — x)^, since the number of sites taking 
part in intcrorbital exchange interaction is reduced by 
a factor of 1 — a; . In the materials with Jahn- Teller 
ions, the orbital Hamiltonians are more complicated than 
the Heisenbcrg-likc Hamiltonian considered in this paper. 
Preliminary calculations for the Hamiltonian correspond- 
ing to the superexchange mechanism of the orbital order- 
ing^ show that the obtained results remain qualitatively 
the same. However, in real substances, there also exists 
a possibility of 00 without local distortions, correspond- 
ing to complex combinations of eg orbitals^'*, which needs 
a special analysis. 
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